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ABSTRACT 

This  paper,  proposes  a  non-linear  constrained  mathematical  model  to  search  the  optimal  cutting  parameters  such 
as  cutting  speed,  feed  rate  and  depth  of  cut.  The  model  deal  with  the  multi-pass  turning  operation  which  is  comprises  of 
multi-pass  rough  machining  and  finish  machining.  An  optimization  technique  based  on  bio-geography  based  optimization 
has  been  introduced  to  optimize  the  multi-pass  roughing  and  single-pass  finishing  parameters  to  achieve  minimum 
production  cost.  BBO  is  very  effective  and  efficient  method  for  optimization  problem.  An  example  is  adopted  from  the 
literature  to  solve  the  proposed  optimization  problem. 

KEYWORDS:  Multi-Pass  Turning  Operation,  Cutting  Parameters,  Bio-Geography  Base  Optimization 

INTRODUCTION 
Nomenclature 

UC  Unit  production  cost  except  material  cost  ($/piece) 

Cm   Cutting  cost  by  actual  time  in  cut  ($/piece) 

Ci    Machine  idle  cost  due  to  loading  and  unloading  operations  and  idle  motion  time  ($/piece) 
Cr    Tool  replacement  cost  ($/piece) 
Ct    Tool  cost  ($/piece) 

Vr  Vs     Cutting  speed  in  rough  and  finish  machining  (m/min) 

VrL  VrU  Lower  bound  and  upper  bound  of  cutting  speed  in  rough  machining  (m/min) 
Vsl ,VsU  Lower  bound  and  upper  bound  of  cutting  speed  in  finish  machining  (m/min) 
/r,/s  Feed  rates  in  rough  and  finish  machining  (mm/rev) 

/ rL,/ ru  Lower  bound  and  upper  bound  feed  rates  in  rough  machining  (mm/rev) 
/sl,/su  Lower  bound  and  upper  bound  feed  rates  in  finish  machining  (mm/rev) 
dr,ds  Depth  of  cut  for  each  passes  of  rough  and  finish  machining  (mm) 

dlL,drU  Lower  bound  and  upper  bound  depth  of  cut  for  each  passes  of  rough  machining  (mm) 
dsL,dsU    Lower  bound  and  upper  bound  depth  of  cut  for  each  passes  of  finish  machining  (mm) 
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n     Number  of  rough  cuts  (an  integer) 
dt     Depth  of  material  to  be  removed  (mm) 
D,L  Diameter  and  length  of  work  -piece  (mm) 
kQ    Direct  labor  cost  +  overhead  ($/min) 
kt     Cutting  edge  cost  ($/edge) 

tmr,tms,tm  Rough,  finish  and  actual  machining  times  (min) 

te,tr  Tool  exchange  and  tool  replacement  times,  (min) 

h!,h2  Constants  related  to  tool  travel  and  approach  /  departure  time  (min) 

T,Tr,Ts  Tool  life,  expected  tool  life  for  rough  machining,  expected  tool  life  for  finish  machining  (min) 

Tp    Tool  life  of  weighted  combination  of  Tr  and  Ts  (min) 

TlT\j     upper  and  lower  bounds  for  tool  life  (min) 

p,q,r,C0  Constants  of  tool  life  equation 

Fr,Fs  Cutting  forces  during  rough  and  finish  machining  (Kgf) 

ki,(a,v    Constants  of  cutting  force  equation 

P,,PS  Cutting  power  during  rough  and  finish  machining  (KW) 

Pu    Maximum  allowable  cutting  power  (KW) 

SC   Limit  of  stable  cutting  region  constraint 

Qr,Qs     Chip  tool  interface  rough  and  finish  machining  temperatures  (C°) 
Qu   Maximum  allowable  chip  tool  interface  temperature  (C°) 
OBJECTIVE 

Optimization  has  become  very  important  in  this  competitive  world.  In  any  field  we  can  observe  that  optimization 
has  become  a  key  point  to  success.  In  this  paper,  optimization  of  turning  operation  is  done.  Turning  process  is  the  most 
common  machining  process  for  cutting  operation.  The  process  of  metal  removal  using  turning  operation  comprise  of  two 
stages  e.g.  rough  machining  stage  and  finish  machining  stage.  Several  variables  such  as  feed  rate,  cutting  speed,  depth  of 
cut,  work  material,  its  properties  and  characteristics  of  output  variables  such  as  production  cost,  production  time,  no  of 
iteration,  tool  life,  surface  roughness,  temperature,  cutting  force  etc.  should  be  considered  to  get  the  final  products  that 
meet  the  specification.  Many  traditional  mathematical  programming  techniques  have  been  used  to  solve  optimization 
problems,  but  these  techniques  had  so  many  drawbacks.  These  techniques  couldn't  solve  the  multimodal  problems  as  they 
give  only  local  optimal  solutions  and  cannot  solve  the  problem  having  so  many  constraints.  Metaheuristic  techniques  such 
as  particle  swarm  optimization,  genetic  algorithm,  ant  colony  algorithm,  evolutionary  algorithm,  bio-geography  based 
algorithm  are  being  utilized  to  solve  these  optimization  problems.  In  this  paper,  Bio-geography  based  algorithm  is  being 
used  to  optimize  the  cutting  parameters  of  multi-pass  turning  process.  Bio-geography  based  optimization  is  a  nature  based 
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algorithm.  BBO  is  the  study  of  geographical  distribution  of  biological  organism  such  as  animal,  plants  over  time  and  space. 
The  main  objective  of  this  process  planning  is  to  find  the  appropriate  cutting  parameters  which  can  provide  max  profit  to 
the  company  and  could  be  available  to  the  customer  with  appropriate  quality  and  lead  time.  Crookall  and 
Venkataramani[l]  studied  the  optimization  technique  for  the  turning  operation. 

A  probabilistic  approach  has  been  developed  to  the  determination  of  the  optimum  cutting  conditions.  But  in  this 
paper  constraints  are  not  considered.  Then,  Shin.  Y.  C.  and  Joo  Y.S.[3]  Optimize  the  machining  conditions  with  practical 
constraints.  They  proposed  the  mathematical  formula  to  optimize  the  cutting  parameters  by  dynamic  programming 
approach.  Then,  Onwubolu  and  Kumalo  [4]  approach  with  Optimization  of  multipass  turning  operations  by  genetic 
algorithms.  They  optimize  the  same  problem  using  Genetic  algorithm  and  proved  more  efficient  as  compare  to  the 
Simulated  Annealing.  Vijayakumar  [5]  use  Ant  Colony  algorithm  to  optimize  and  proved  their  method  give  better  result 
with  less  iteration  and  within  short  time.  They  compare  their  result  with  other  algorithm,  and  come  to  a  conclusion  that 
ACO  is  very  effective  method.  Wang  YC.  [6]  works  on  optimization  of  Ant  colony  algorithm  and  find  that  paper  of 
vijaakumar[5]  where  he  did  not  prove  the  optimal  values  they  found  for  the  depth  of  the  rough  cuts  and  the  finishing  cut 
and  the  constraint  related  to  the  no  of  cuts.  Abderrahim,  and  imane  [7]optimize  the  cutting  parameters  by  newly  developed 
firefly  optimization  technique.  Tsai  [12]  came  with  simulated  annealing  approach  for  optimization  of  multi-pass  turning 
operations  They  solve  the  optimization  problem  using  same  mathematical  formulae  by  Simulated  Annealing  algorithm. 
This  algorithm  giving  less  production  cost  but  the  time  consumption  is  more  than  other. 

The  algorithm  is  considered  in  the  continuous  constraints  optimization  problem  where  the  task  is  to  minimize  the 
production  cost.  The  researcher  compared  with  other  algorithm  such  as  genetic  algorithm,  simulated  annealing  and  found 
that  the  result  obtain  near  optimal  values.  In  this  paper,  Bio-geography  Based  optimization  technique  (BBO)  is  used  to 
optimize  the  multipass  turning  operation.  We  are  using  the  same  mathematical  formulation  given  by  Shin,  Y,  C,  and  Joo, 
Y.S.(3) 

MULTI-PASS  TURNING  MODEL 

The  main  objective  of  the  paper  is  to  find  the  optimal  cutting  parameter  including  cutting  speed,  feed  rate,  depth 
of  cut,  number  of  rough  cuts  and  no  of  iteration  in  order  to  minimize  the  unit  production  cost  and  simultaneously  achieve  a 
good  surface  finish  considering  certain  constraints  such  as  roughness,  power  etc.  The  mathematical  model  used  by  Chen 
and  Tsai  [12]  is  preferred  adapted  to  find  the  optimal  machining  parameters. 

Cost  Function 

The  unit  production  cost,  UC,  for  the  multi-pass  turning  operations  problem  consist  of  four  components. 

•  Machining  cost  by  actual  time  in  cutting  operation,  Cm. 

•  Machine  idle  cost  due  to  loading  and  unloading  operations  and  idle  tool  motion,  Q 

•  Tool  replacement  cost,  CR. 

•  Tool  cost,  CT 
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Unit  Production  Cost 

Using  equations  (1),  (2),  (3)  and  (4),  the  unit  production  cost,  UC,  is  defined  as: 
UC  =  CM  +  Ci  +  CR  +  CT; 

.  r     nDL      ,dt— ds.    ,       nDL     n      r         .    t     .   .  ,dt— ds  ,    . ... 

UC  =  k0[  (  )  +  +  o  tc  +  (h,L+h2)(  +  1)  +  (5) 

,    ter     nDL      .dt— ds.    ,       nDL     -.      Ktr     nDL      ,dt— ds.    ,       nDL  -. 

kn—[  C  )  "I  ]  +  — [  C  )  "I  ] 

Tp  1000Vr/r  v   dr   J      lOOOFs/s       Tp  1000Vr/r v   dr   J  lOOOFs/s 

Surface  Quality  Function 

The  surface  quality  function  is  described  as: 

Rakvsxlfsx2dsx3  (6) 

Where  x1,x2,x3and  k  are  the  constants  relevant  to  specific  tool-  work  piece  combination.  Since  reaching  an  idle 
surface  finish  is  not  necessarily  required  in  industry,  so  in  this  paper  surface  quality  function  is  treated  as  a  constraint  and 
therefore  constraint  satisfaction  method  is  adopted 

Cutting  Condition  Constraints 

Several  cutting  constraints  are  considered  during  roughing  and  finishing  operations  as  follows: 
Parameter  bounds 
Tool-life  constraint 
Cutting  force  constraint 
Power  constraint 
Surface  finish  constraint 
Stable  cutting  region  constraint 
Chip-tool  interface  temperature  constraint 
Surface  finish  constraint  (for  finishing  stage  only) 
Roughing  and  finishing  parameter  relations 
Depth  of  cut  equality  constraint. 
Roughening  Machining 
•     Parameter  Bounds 

Due  to  the  limitations  on  the  machine  and  cutting  tool  and  due  to  safety  of  machining  the  cutting  parameters  are 
limited  with  the  bottom  and  top  permissible  limit: 

Cutting  speed:  VrI<Vr<Vru' 

Feed  rate:  /rJ</r</rU  (7) 
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Depth  of  cut:  drI<dr<dru 

•  Tool-Life  Constraints 

The  constraint  on  the  tool  life  is  taken  as 

Ti<rr<ru  (8) 

•  Cutting  Force  Constraint 

The  maximum  amount  of  cutting  forces  should  not  exceed  a  certain  value  as  higher  forces  produce  shakes  and 
vibration. This  constraint  is  given  by 

F1-k,(/,r(ar)°<Fu  (9) 

•  Power  Constraint 

The  nominal  power  of  the  machine  Pv  limits  the  cutting  process 

FrVr 
61201] 

Efficiency  q=0.85 

•  Stable  Cutting  Region  Constraint 

This  constraint  is  given  by 

(V/frXd^SC  (11) 

•  Chip  Tool  Interface  Temperature  Constraint 

The  temperature  generated  during  the  cutting  operation  should  not  exceed  the  permissible  limit.  It  is  expressed  as; 

Qr=ka(V1)T(fr),(4)*<Qu  (12) 
Finish  Machining 

All  the  constraint  other  then  the  surface  finish  constraint  are  similar  for  rough  and  finish  machining 

•  Surface  Finish  Constraint 

The  surface  finish  constraints  is  given  by: 

£<(SR)u  (21) 
The  cutting  parameter  relation  constraints  are: 

Vs>k3Vr  (22) 
f^k,fs  (23) 
d^k5ds  (24) 
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•     The  Number  of  Rough  Cuts 

The  possible  number  of  rough  cuts  is 

dt-ds 

n=—  (25) 
The  restriction  is: 

dt-dsl      dt-ds  dt-dsU 

 <   <    (26) 

drl  dr  drU 

The  optimization  problem  in  multipass  turnings  can  be  divided  into  m=(nu-nL+l)  sub  problems,  in  each  of  which 
the  number  of  rough  cuts  n  is  fixed.  So  the  solution  of  the  whole  optimization  problem  is  divided  into  searching  the 
optimal  result  of  m  sub  problems  and  the  minimum  of  them  is  the  objective  of  the  whole  optimization  problem. 

BIO-GEOGRAPHY  BASED  OPTIMIZATION 

Bio-geography  based  optimization  [10]  is  the  advance  of  evolutionary  algorithm  used  for  searching  global  optima. 
It  is  inspired  by  nature  and  is  followed  by  immigration  and  emigration  of  species  from  one  habitat  to  the  other  in  search  of 
rich  habitats.  The  solution  is  represented  by  habitat  with  Habitat  Suitability  Index.  This  process  is  consists  of  two  process 
migration  and  mutation. 

In  migration  process,  many  parents  contribute  to  a  single  offspring.  Migration  is  used  to  change  the  existing 
solution  and  modify  the  solution.  Migration  is  the  probabilistic  process  that  adjusts  the  habitats  X;.  The  probability  is 
modified  w.r.t  immigration  rate  and  the  source  of  modified  probability  comes  from  Xj  proportional  to  the  emigration  rate 
and  is  expressed  as: 

^=1(1-—)  (27) 

v    Smax  '  ' 

U=—  (28) 
r   Smax  ' 

Where,  I  is  the  maximum  possible  immigration  rate,  E  is  maximum  emigration  rate,  s  is  no  of  species,  Smax  is 
max  no  of  speci 

Mutation  is  the  second  stage  of  bio-geographical  process.  During  mutation  drastic  change  made  to  the  HSI  of 
habitat  due  to  certain  events.  Mutation  increases  with  diversity  among  the  candidate  size.  Candidate  solutions  is  depend  on 
the  mutation  probability  as  given  in  below. 

M(s)=(l-Ps)Mmax/Pmax  (29) 

Where  Mmax  is  a  user-defined  parameters,  Ps  is  the  species  count  of  the  habitat,  Pmax  is  maximum  species  count. 
Mutation  is  done  on  the  basis  of  mutation  probability  of  each  habitats  and  by  replacing  the  existing  suitability  index 
variable  (SIV)  of  habitat  with  the  another  generated  suitability  index  variable. 
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Choose  parameter  n,  iteration  P„ 


Randomly  generated  SI  V  of  each  island  that 
form  island 

 1  


J  =  1  +  1 


Evaluate  SI  for  each  island  and  identify  the  island  having  best  fitness  value 


Perform  migration  probabilistically  on  those  SIV 


Perform  Mutation  probabilistically'  on  those  SIV  based  on  p„ 


SIV  of  b  est  island  is  optimal  solution  and  apply  on 
solution 


start 


Figure  1:  Flow  Chart  of  BBO 


RESULTS 


An  example  is  considered  from  [12]  to  demonstrate  and  validate  optimization  problem  though  turning  process.  In 
this  section,  a  numerical  study  is  performed  to  analyze  the  effects  of  decision  variables.  Evaluation  and  comparison  is  also 
done  with  the  results  given  in  the  previous  literature 

Table  1:  Machining  Data 


Parameter 
/Constrain 
ts 

Value 

Parameter/ 
Constraint 

Value 

Parameter/ 
Constraints 

Value 

Parameter/ 
Constraints 

Value 

D 

50  mm 

L 

300  mm 

Co 

6*10" 

d, 

6mm 

vRl, 

500  m/min 

vrI 

50  m/min 

vsU 

500  m/mn 

vsL 

50  m/min 

frU 

0.9  mm/rev 

0.1  mm/rv 

fsU 

0.9  mm/rv 

fsL 

0.9  mm/rev 

drU 

3.0 

1.0 

dsU 

3.0 

dsL 

1.0 

P 

5 

q 

1.75 

r 

0.75 

ri 

0.85 

K 

0.5  $/min 

k, 

108 

k2 

132 

k3 

1.0 

U 

2.5 

k5 

1.0 

kt 

2.5  $/ege 

e 

0.8 

T 

'  min 

25  min 

T 

A  max 

45  min 

F 

1  max 

200  kgf 

p 

A  max 

5  kw 

0.75 

V 

0.95 

SRu 

10  um 

R 

1.2  mm 

X 

2 

V 

-1 

SC 

140 

Qu 

1000  °c 

X. 

0.4 

0.2 

6 

0.105 

to 

0.75  min/piece 

hi 

7*10-4 

h2 

0.3 

t„ 

0.75  min/piece 

te 

1.5  mm/ege 
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Table  2:  The  Optimized  Turning  Parameters 


Cutting  Parameters 

Cost 

Iteration 

SPEED 

157.3932 

ROUGHGING 

FEED  RATE 

0.1142 

DEPTH  OF  CUT 

2.6482 

SPEED 

458.4504 

FINISING 

FEED  RATE 

0.8641 

DEPTH  OF  CUT 

2.8079 

1.2086 

42 

Table  3:  Results  of  Optimization  Using  Different  Algorithm 


Algorithm 

Cost 

BBO 

1.2086 

ACO 

1.8450 

BFOA 

1.4231 

FIREFLY 

2.5271 

DE 

2.8680 

SA 

2.2959 

GA 

1.761 

Graph 


Figure  4:  Graph  of  Minimum  Cost  Obtained  Vs  Iteration 

Here  in  this  problem,  certain  assumptions  are  taken  as  follows: 
Population  size  =  100; 
Probability  of  modification  =  1 ; 
Probability  of  mutation  =  0.3; 
Maximum  mutate  =  0.6; 

Table2  Showing  the  optimal  cutting  parameters  of  the  turning  process  used  to  get  minimum  production  cost.  From 
the  table  we  can  observe  that  bio-geography  based  optimization  is  giving  minimum  production  cost  e.g  1.2086  which  is 
minimum  than  the  algorithm  proposed  by  various  researchers.  Ant  colony  algorithm  (ACO)  proposed  by  vijaykumar  [9] 
achieves  production  cost  equal  to  1.8450  and  also  violates  so  many  constraints.  While  simulated  annealing  proposed  by 
Onwubolu  [11]  arrives  at  a  production  cost  (UC)  value  of  1.761.  In  the  above  table  we  can  observe  that  the  proposed 
algorithm  insures  a  minimum  production  cost  with  45  %  reduction  in  the  production  cost  then  in  firefly  algorithm.  Whereas 
in  other  algorithm  the  production  cost  is  higher  as  compare  to  the  proposed  algorithm.  In  GA  algorithm,  the  production 
cost  is  less,  but  time  taken  during  iteration  is  higher  as  compare  to  proposed  algorithm.  Hence,  the  proposed  BBO 
algorithm  is  effective  to  solve  such  kind  of  optimization  problem 
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CONCLUSIONS 

In  this  paper,  problem  related  to  optimization  of  cutting  parameters  of  turning  process  considering  constraint  is 
implemented  and  solved  by  using  bio-geography  based  algorithm.  The  main  aim  of  this  project  is  to  minimized  total 
production  cost.  The  cutting  parameters  as  cutting  speed,  feed  rate,  depth  of  cut  during  roughing  and  finishing  pass  are  the 
main  process  parameters  whose  optimal  values  effects  the  machining  process.  Many  researchers  solve  this  problem  taking 
various  algorithms,  but  total  production  cost  (UC)  came  in  this  paper  is  comparatively  low  as  compare  to  other  algorithm 
implemented  by  other  researchers.  BBO  provides  better  result  with  quality  and  feasibility.  So  we  can  say  BBO  is  a  reliable 
and  efficient  method  for  solving  complex  machining  program 
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